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THE  AERODYNAMIC  HEATING  OF  A  COMPOSITE  FLAT  PLATE* 


Prepared  by: 
James  L.  Rand 


ABSTRACT:  The  purpose  of  this  study  was  to  find  solutions  to 

the  heat-conduction  equation  in  both  one  and  two  dimensions 
when  applied  to  a  composite  flat  plate.  The  two-dimensional 
model  is  that  of  two  slabs  with  different  thermal  properties, 
in  good  thermal  contact  with  each  other  while  being  exposed  to 
the  same  boundary  layer.  The  solution  takes  the  form  of 
Fourier  series  which  has  been  found  to  be  orthogonal  with 
respect  to  a  "weighting"  function  identical  to  the  "weighting" 
function  found  by  Mayer  in  the  one-dimensional  analysis  of  a 
composite  slab.  The  Van  Driest  turbulent  boundary-layer 
theory  has  been  applied  to  the  one-dimensional  problem  to 
determine  the  constant  film  coefficient  boundary  condition. 

It  has  been  found  that  for  those  values  of  the  Fourier  modulus 
less  than  0.1,  a  maximum  error  in  surface  temperature  ratio  of 
8  percent  may  be  expected  while  the  temperature  distribution 
in  general  is  unaffected  for  the  case  considered. 
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INTRODUCTION 

In  recent  years  attention  has  been  focused  on  the  problems 
associated  with  the  extremely  high  temperatures  encountered  by 
a  vehicle  during  re-entry.  A  working  knowledge  of  the  nature 
of  the  heat -convection  problem  Is  required  by  the  engineer  so 
that  the  vehicle  may  be  designed  to  function  properly  In  such 
a  hlgh-temperature  environment.  As  a  result,  the  rate  of  heat 
transfer  to  the  vehicle  becomes  an  extremely  Important  design 
parameter.  However,  since  the  nature  of  all  aerodynamic 
heating  problems  Is  not  yet  fully  understood,  most  engineers 
must  resort  to  experimental  data  for  the  prediction  of  the 
various  heating  parameters. 

Experimental  data  may  be  obtained  by  various  methods,  such 
as,  wind  tunnel  or  shocktube  tests.  In  both  types  of  testing, 
an  attempt  is  made  to  simulate  both  the  geometry  and  environment 
the  full  scale  vehicle  might  be  expected  to  experience.  The 
heat -transfer  parameters  are  then  calculated  from  some  knowledge 
of  the  temperature  history  of  the  model. 

One  of  the  more  popular  models  (fig.  l)  used  to  obtain 
experimental  data  is  the  smooth  flat  plate  placed  in  the 
direction  of  a  high-temperature  uniform  flow  in  which  the 
pressure  is  constant.  The  temperature  distribution  in  the 
plate  is  then  given  by  the  theory  of  heat  conduction,  where  the 
boundary  conditions  at  the  surface  exposed  to  the  flow  are 
given  by  the  various  boundary-layer  theories.  Since  this 
provides  one  of  the  simplest  of  mathematical  models  it  is 
possible  to  study  the  nature  of  the  phenomenon  more  accurately. 

Normally,  the  heat-conduction  problem  may  be  considered  to 
be  one  dimensional  with  boundary  conditions  of  varying  degrees 
of  complexity.  For  relatively  low  rates  of  heat  transfer  or 
short  periods  of  heating  (ref.  (1)),  the  flat  plate  may  be 
considered  to  be  a  semi-infinite  solid.  This  arises  from  the 
fact  that  the  total  heat  flux  is  not  sufficient  to  penetrate 
the  plate  during  the  period  of  the  test.  As  the  thickness  of 
the  plate  is  reduced,  the  mathematical  model  must  be  changed  to 
that  of  a  finite  slab  with  various  boundary  conditions  on  the 
rear  surfaces.  As  a  practical  application,  the  temperature 
distribution  in  a  calorimeter  gage  (ref.  (2))  is  obtained  when 
this  boundary  condition  is  that  of  zero  heat  flux.  When  the 
thickness  of  the  plate  is  reduced  still  further,  the  temperature 
may  be  considered  constant.  Utilizing  this  boundary  condition, 
the  response  of  a  thin  film-resistance  thermometer  may  be 
obtained.  Carslaw  and  Jaeger  (ref.  (3))  have  described  most  of 
these  conditions  and  their  corresponding  solutions  by  various 
methods  in  their  text. 
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The  shocktube  Is  often  used  to  produce  extremely  high 
rates  of  heat  transfer  because  of  Its  ability  to  generate 
extremely  high  enthalpies  and  temperatures.  Heat -transfer 
rates  may  be  obtained,  as  usual,  by  Imbedding  a  calorimeter 
gage  In  the  flat  plate  exposed  to  the  flow.  However,  due  to 
the  extremely  high  pressures  associated  with  this  type  of  test, 
severe  problems  are  encountered  In  bonding  and  sealing  the  gage 
Into  position.  As  a  result,  the  gage  must  be  made  much  thicker 
than  normal  and  the  resulting  area  In  thermal  contact  with  the 
bonding  material  becomes  quite  large.  The  heat  flow  in  the 
gage  may  become  two  dimensional  depending  on  the  thermal 
properties  of  the  bonding  material.  In  order  to  determine  the 
temperature  distribution  in  the  gage  the  flat  plate  is  assumed 
to  be  formed  by  two  dissimilar  materials  in  good  thermal  contact 
as  shown  in  figure  1.  Since  both  materials  are  exposed  to  the 
same  boundary  layer  and  to  each  other,  a  solution  to  the  heat- 
conduction  equation  in  two  dimensions  must  be  developed.  With 
this  solution,  the  effect  of  the  aerodynamic  heating  of  a 
composite  flat  plate  may  be  calculated. 


ONE-DIMENSIONAL  ANALYSIS 

It  has  been  mentioned  in  the  preceding  section  that  the 
temperature  distribution  of  a  calorimeter  gage  may  be  calculated 
from  the  theory  of  heat  conduction  in  a  slab  of  finite 
thickness.  If  the  flow  of  heat  may  be  considered  one  dimen¬ 
sional,  the  mathematical  model  used  for  this  calculation  is 
shown  in  the  following  sketch. 

L-L-L_i— L-i-l-l— i_i— 1— i  — 

y=Q 
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The  differential  equation  governing  the  flow  of  heat  Is  given 
by: 

=  0  (1) 

Where  T  =  T  (y,t)  Is  the  temperature  above  the  Initial 
temperature  of  the  gage  (i.e.  T  =  0  when  t  =  o).  At  times 
greater  than  zero,  the  temperature  at  the  surface  Is  governed 
by  the  equation: 

k£|  +  h[T.w  -TU.t)]  =  0  i  >0  (a) 

°Y  ly-'jt  1 

The  flux  is,  therefore,  proportional  to  the  difference  between 
the  surface  temperature  and  some  temperature  determined  by  the 
surrounding  medium.  The  proportionality  constant  "h"  is  the 
film  coefficient  of  heat  transfer  which  Is  assumed  constant  for 
this  analysis. 

At  the  back  surface  the  heat  flux  is  assumed  to  be  zero, 
therefore: 


in./0  (3> 

Let  us  assume  a  solution  of  the  form: 

T(y.t)  =g(y)e',,t  +  C  <*) 

Upon  substitution  into  equation  (1)  we  obtain  a  differential 
equation  in  g(y)  which  has  as  its  solution: 

g(yj  =  A  cos  +  Bsm  J =r  (5) 

By  combining  equations  (^)  and  (5)  and  applying  the  boundary 
condition  at  y  =  0,  we  see  that: 


6  =  0  (6) 
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Upon  applying  the  boundary  condition  at  y  =  1  we  have,  after 
some  rearrangement : 


0  tAN  '  L  *  0 


(7) 


and 


C  —  J  AW 


where  fi- 


Xi 

f>r 


and  L  = 


u_ 

K  ' 


the  Biot  modulus 


Then 


T(y.t)=T^ 


cos 


(8) 


(9) 


where  the  A'8  are  the  positive  roots  of  equation  (7),  and 
f *  ,  the  Fourier  modulus,  a  non-dimensional  time. 

The  coefficients  An  must  be  determined  from  a  Fourier  analysis 
of  the  series  at  time  zero.  Therefore: 


»  S  | 


A*  cos 


By  Fourier  analysis 


or  finally: 


Ay 

i 

Ta*j 

t-0 

(10) 

1: 

Jt 

=  -X  , 

loos  4-V, 

(11) 

\coS*  A  1  Jy 
»  X  7 

s  -2  Xiu* 

(12) 

This  solution  is  a  special  case  of  the  solution  obtained  by 
Carslaw  and  Jaeger  to  the  problem  of  the  finite  slab  with 
arbitrary -boundary  conditions. 


A  more  exact  analysis  of  this  problem  is  sometimes 
achieved  by  taking  into  account  the  loss  of  heat  at  the  back 
surface.  This  is  done  by  assuming  a  second  dissimilar  material 
to  be  in  good  thermal  contact  with  the  back  of  the  gage.  At 
the  interface  of  the  two  materials  both  the  temperature  and 
heat  flux  must  be  identical,  therefore: 


T, 


(13) 
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•“»  K,  •§£  -  K,  ^7  U*> 

The  second  material  may  then  be  thought  of  as  a  semi-infinite 
solid  or  another  finite  slab  as  shown  in  the  following  sketch. 

..  i  t  <  i  <<♦  +  <*<<  < _ 

y=*2 


y*li 


y  =  o 

-00  - \\\\\\\xvx\  \\\\ - +0° 

This  problem  has  been  treated  by  many  authors  since  it  is  an 
excellent  mathematical  model  for  a  variety  of  engineering 
problems.  In  most  treatments  of  this  problem  the  film 
coefficient  is  again  considered  to  be  constant.  Mayer 
(ref.  (4))  has  obtained  a  solution  by  the  methods  of  Fourier 
analysis  which  is  very  useful  for  large  values  of  time. 
Wasserman  (ref.  (5))  and  Campbell  (ref.  (6))  have  obtained 
solutions  for  small  values  of  time  for  this  problem. 

If  the  heating  rate  or  testing  time  is  extremely  short, 
the  surface  temperature  will  not  change  appreciably  and  the 
boundary  condition  given  by  equation  (2)  may  be  approximated 
by: 

-x  17U/  Q‘  (15) 

where  Q0  is  the  heat  flux  at  time  zero,  and  is  considered  to 
be  constant.  Rose  (ref.  (2))  has  used  this  boundary  condition 
and  solved  the  problem  by  LaPlace  transform  methods  to  obtain 
short-time  solutions  which  normally  do  not  have  the  convergence 
problems  associated  with  the  Fourier  analysis. 
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TWO-DIMENSIONAL  ANALYSIS 


The  one-dimensional  solutions  previously  discussed  are 
excellent  treatments  of  the  heat- conduction  problem  when  the 
y-dlmension  is  much  smaller  than  the  x-dimension.  However, 
when  these  dimensions  are  of  the  same  order  of  magnitude,  care 
must  be  taken  to  consider  any  possible  end  effects. 

Therefore,  consider  the  composite  flat-plate  system  Shown 
in  figure  1.  If  we  restrict  our  problem  to  that  of  finding  the 
temperature  distribution  for  small  values  of  time,  we  may 
oonsider  at  the  back  surface,  the  heat  flux  in  the  y-direction 
to  be  everywhere  zero  and  at  the  front  surface  the  film  coeffi¬ 
cient  to  be  constant.  Then  by  considering  the  symmetry  of  the 
system, the  heat  flux  at  the  center  of  each  element  is  zero,  and 
we  may  reduce  our  problem  to  that  shown  in  the  following  sketoh. 


y 


The  differential  equations  governing  the  flow  of  heat  are  given 
by: 


and 


-+  xt 


VTt 

V 


0 


(16a) 

(16b) 
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where  T"="T(*.y,t)ls  the  temperature  above  ambient  and  Is  Initially 
zero  throughout.*  The  boundary  conditions  for  this  problem  are 
as  follows: 

In  region  one: 


ill  =  0 

(17) 

&i  = 0 

(18) 

*  Y  1  (ir.  0,  t) 

(19) 

In  region  two: 

151  =  0 

(20) 

■>yl 

(21) 

K  AS 
**  'iy 

+  .<)  -TaJ= 0 

(Xx.l.t)  1  1 

(22) 

And  at  the  interface: 

T  (<*.y 

',+)  =  Tx  Clo.y,*) 

(23) 

1/ hi  --k  u»i 

'  2  }*,  l(b.y.t) 

(24) 

Seeking  solutions  by  separation  of  variables  let  us  assume 
solutions  of  the  form: 


\\ 

T(x.,y,t)  e  +  C, 


(25a) 


♦The  coordinate  system  has  been  selected  as  shown  in  figure  8, 
such  that  at  ^*0,  the  heat  flux  is  zero. 
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and 


Tt  (x»,  y,*)  =  f,  <x>  (y) 


(25b) 


Substitution  of  these  functions  Into  equations  (16)  yield: 


(26a) 

(26b) 


Since  f ( x )  and  g(y)  are  functions  of  different  independent 
variables,  each  part  of  equations  (26)  must  equal  a  constant; 
therefore: 


and  similarly  say, 

then  we  obtain: 


and 

and 


(27a) 

(27b) 


f,U)  -  A,  cos  v-|-  +  B.  sim  v-J- 


(28) 


and  gi(y)=  E,  cos  +  D,  cos  j[  (29) 

and  similar  solutions  for  and  y)  .  It  is 

apparent  from  equations  (17)  and  (18)  that  132=0}=  0. 
Likewise,  it  may  be  seen  that  Bo  =  Do  =  0.  Prom  equations  (19) 
and  (22),  we  may  conclude  that  61=62=  Taw,  so  that  the  time 
dependent  term  can  be  canceled  yielding: 


K,'3fU+h’'<e)  *  0  <«•) 

and  K«  =  0  (22'> 

Substituting  equations  (29)  into  (19')  and  (22')  we  obtain: 


E,  £  -  L,  cos  £ ]  —  0 


(30a) 
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and  [  a  sin  «  -  L2  cos  <*  ]  -  0 

(30b) 

.  M 

Where  L„*  18  again  the  Blot  modulus,  the 
Ei  and  Eg  cannot  equal  zero  lest  the  solution 
dimensional,  therefore: 

constants 
become  one- 

0  t*N  /8  -  L ,  •  0 

(31a) 

and  a  t*  *>  or  -  L2  =  0 

(31b) 

Let  us  now  rewrite  equations  (25)  with  the  help  of  (28)  and 
(29): 

v2 

T(X,y,t)  =A  cos  V-5-  cos^  -j-  c.  +  Taw 

(32a) 

*T^(x,y,t)  =  B  cos  S  -j”  cos  o(-~  c.  ^  +  Tiu 

(32b) 

Prom  equations  (23)  and  (24)  we  have: 

A  c .os  x  =  B  cos  $  cos  «  j-  (23') 

“A  K,  sim  cos  /&  -j-  -  8  K2  i  cosof-i.  (24') 


Dividing  (24')  by  (23')  we  obtain: 

K,  b  y  tAiay  =  -  a  S  tAw  &  (33) 

Prom  equations  (26)  and  (27)  we  may  write: 

*•  [(i)*+  (ft  ]  *  *.  [ff)!+  fi-f]  <3“> 

Ve  may  now  write  our  solution  In  Its  final  form: 
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and 


T  (x.y. 0  *  Taw  ♦|EAM(I  cos  Co* 

y, t)  =  "Taw  Bm-  cosSMM-r'  casctAe‘ *““* 

#*fcl  N*»  D  X 


(35a) 

(35b) 


where  the  aw»  and&,'s  are  the  positive  roots  of  equations  (31)# 
and  the  YmN‘s  and  &  !•.„'&  are  determined  from  the  simultaneous 
solution  of  equations  (33)  and  (3^).  The  AlMU’s  may  be 
determined  from  (26a)  and  (27a)  as: 


x*»  -  tI(4^1*^2] 


(36) 


The  coefficients  of  the  series  indicated  in  equations  (35) 
must  be  determined  from  a  Fourier  analysis  of  the  entire 
composite  flat-plate  system  at  time  zero. 

In  order  to  expand  the  initial  temperature  into  a  Fourier 
series,  some  orthogonality  relation  must  be  determined.  It 
has  been  shown  in  Appendix  A  that: 

f,c,[eos  V^-^cos  A*j][cos  VP1-J  cos /3p  j] due//  + 

► 

l£ek[tos*«-'b  c°s«hx][C0*  Inf co*  dy  =  0  (37) 


0  0 


Therefore,  evaluating  equation  (35)  at  time  zero,  multiplying 
by  the  orthogonal  function,  and  integrating  over  the  entire 
composite-wall  system  we  obtain: 

JV,  j  t 

ft  c,[cos  cos  dx 

•  d  J 

«  (JC.fcosY.Jf  dy  ^tf:»iLfcrt9m^(Jr(38) 


From  equation  (23')  we  may  write: 

COS  Yw- 

•  - 7 - y  A 

COS  O  ns  COX  0(_  — 


b 


(39) 
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Substituting  equation  (39)  Into  (38)  yields  the  following: 

<co* 


a-.[p.  *•  c,~"  *  fciiir)  ftc*  c,~l = 


iC os 

-Ta-[,o  c  D,Mm  + 


t  * 


OP 


•  c 

M  4 

ft 

•• 

4> 


(mM  r  c  c 

(40) 

,-~a  CoS  ft**’]  4*^1 

(41) 

«--b  CftS 

(42) 

(43) 

«.“b  to8  «l*dy 

(44) 

Now  performing  the  Integration  we  have: 


(■'•mm  ~ 


* 


Dimw  “ 


aJ 

1 

<U.  ] 

aJ2 

Y*.  ’ 

■3TC 

SI**  ' 

& m  im* 

Km  CIS  K 


~jj 


(41  ■) 
(42*) 
(43') 
(44‘ ) 


Finally,  the  temperature  Is  obtained  by  substituting  equation 
(36)  Into  (35)  to  obtain: 

T(*,*t)=TA»«f  f  K.c»  ^  *■  (35*.) 

T  X  t 

where  t,*  ~JX  »  18  again  the  Pourier  modulus. 

The  coefficients  n  are  obtained  from  equation  (40)  after 
some  rearrangement  as:  * 
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Amw  l  ^ 


AN* 


4  sin  /$„ 


siw  + 


<J*s 

(A/  C,  6^  J  _ 


nr„„  *  s"i~e*,0 


i.45) 


One  possible  method  of  determining  the  validity  of  this 
solution  is  to  allow  the  material  properties  in  region  two  to 
approach  the  material  properties  in  region  one.  When  this  is 
done  the  composite  flat  plate  behaves  as  a  single  slab  and  the 
above  solution  should  reduce  to  the  one-dimensional  solution 
previously  obtained. 


It  can  be  seen  from  equations  (31)  that  if  Lj  =  I<2»  then 
^n  =^m.  Applying  this  condition  to  equation  (34*  see  that 
by  =*aS.  However,  since  a  and  b  may  take  on  any  values, 

[  »  l«  0.  Substituting  these  values  into  equations  (35'a)  and 
(45)  we  obtain: 


and 


T  ( X,  y,t)  =Ta*  +,£  A  „„  cos^  7-  e 
A  _  _  T  r  4  S 1 0  0m 


y 


AW 


/3~  +  &OS  0„ 


(46) 

(^7) 


which  is  in  fact  the  one-dimensional  solution  previously 
obtained. 


DETERMINATION  OP  EIGENVALUES 

The  temperature  distribution  in  a  composite  flat  plate  le 
a  function  of  x,  y,  and  t  and  is  obtained  by  evaluating  the 
double  Fourier  series  derived  in  the  preceeding  section.  The 
indexing  scheme  is  based  on  the  real,  positive  roots  of  the 
transcendental  equations  (31).  These  equations  have  been 
studied  in  great  detail  by  a  number  of  authors.  Carslaw  and 
Jaeger  (ref.  (3)),  for  example,  have  tabulated  the  first 
6  roots  of  these  equations  for  a  wide  range  of  values  of  the 
Blot  modulus,  "L". 

However,  in  order  to  evaluate  equations  (35'a),  a  third 
index,  1,  must  be  employed  since  the  simultaneous  solution  of 
equations  (33)  and  (34)  does  not  yield  unique  values  for  Y**.j 
or  For  the  purpose  of  illustration,  these  equations  have 

been  presented  in  figure  2,  assuming  the  physical  constants  of 
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the  composite  syBtem  to  be  known.  Using  only  the  first  root  of 
(31a)  and  the  first  6  roots  of  (31b),  equation  (34)  has  been 
evaluated  and  Is  presented  as  the  family  of  curves  represented 
by  solid  lines.  Equation  (33)  Is  presented  as  the  family  of 
curves  represented  by  the  broken  lines.  The  Intersection  of 
these  two  families  represents  the  simultaneous  solution  of  these 
equations.  It  is  obvious  that  any  number  of  values  may  be 
obtained  for  and  which  will  satisfy  these  equations. 
Therefore,  it  is  necessary,  when  evaluating  equations  (35a)  and 
(45),  to  introduce  a  third  indexing  scheme  in  order  to  distinguish 
between  the  intersections  of  each  of  the  solid  curves  shown  in 
figure  2. 


HEAT  TRANSFER  THROUGH  THE  BOUNDARY  LAYER 

As  has  been  previously  mentioned,  the  temperature 
distribution  through  a  flat  plate  is  dependent  primarily  on  the 
film  coefficient  of  heat  transfer.  In  fact,  it  is  common 
practice  to  deduce  this  coefficient  by  measurement  of  the 
temperature  distribution;  therefore,  some  information  as  to  the 
character  of  this  parameter  should  be  reviewed. 

Heat  transfer  to  the  flat  plate  is  caused  by  the  presence 
of  a  boundary  layer  which  was  first  reported  by  Prandtl  (ref.  (7)) 
in  1904  and  today  has  become  the  most  important  branch  of  the 
fluid  dynamics  of  viscous  flow  (ref.  (8)).  Many  researchers 
have  considered  this  phenomenon;  in  fact,  Eckert  (ref.  (9)) 
cites  439  separate  contributions  to  the  heating  problem  alone. 

We  will  restrict  our  consideration  to  the  "turbulent"  boundary 
layer  since  this  produces  the  most  severe-heating  conditions 
during  high-speed  flight.  The  presence  of  this  boundary  layer 
produces  a  shearing  force  at  the  surface  of  the  flat  plate 
given  by: 


_ 

~~  \  y 


(48) 


where  is  the  coefficient  of  visccsity  at  the  surface. 

A  coefficient  of  friction  is  then  defined  on  the  basis  of 
this  quantity  and  the  free-stream  conditions  as: 

2TW 

~  P  uJ  (49) 
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The  heat  flux  at  the  surface  Is  given  by  the  theory  of 
heat  conduction  as: 


w 


(50) 


Finally  the  film  coefficient  Is  defined  In  terms  of  the 
heat  flux  as: 


h  ^ 


^  |j 

where  "Ta«  =T»(  1  +  ^1  and  i\  Is  a  recovery  factor.  This  Is 

essentially  the  temperature  the  surface  might  be  expected  to 
attain  when  Insulated. 


CG--Z) 


(51) 


It  haB  been  reported  (ref.  (10))  that  the  Van  Driest  tur¬ 
bulent  boundary -layer  theory  (ref.  ( 11 ) )  for  a  flat  plate 
predlots  local  wall  friction  coefficients  which  agree  well  with 
experimental  data.  Therefore,  this  theory  will  be  used  as  an 
Illustration  of  the  dependence  of  the  film  coefficient  on 
various  aerodynamic  parameters.  Assuming  a  Prandtl  number  of 
unity.  Van  Driest  has  derived,  along  with  others,  that  the  heat 
transfer  at  the  wall  may  be  related  to  the  shearing  force  as 
follows: 


h  =  (52) 

Then  based  on  Prandtl* s  mixing-length  theory  and  knowledge  of 
certain  limiting  conditions,  Van  Driest  obtained  his  final 
equation  for  heat  transfer  through  a  turbulent- boundary  layer 
in  terms  of  the  local  coefficient  of  friction: 


0*1  -*•  ,09,o  R~c*.  "  0.76  los,. 

where  A*=  J  and  B-^l+TfO^E  "  i 

Therefore,  the  above  equation  allows  the  coefficient  of  friction 
and  therefore  the  film  coefficient  of  heat  transfer  to  be 
computed  for  any  arbitrary  wall- temperature  ratio,  Reynolds 
number  and  free-stream  Mach  number. 
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It  should  be  noted  that  the  analyses  In  the  preceding 
sections  were  based  on  the  assumption  that  the  film  coefficient 
of  heat  transfer  remains  constant  during  the  period  of 
Investigation.  However,  It  has  been  shown  that  this  coefficient 
is  a  function  of  not  only  the  aerodynamic  parameters  of  Mach 
number  and  Reynolds  number  buv  also  the  surface  temperature  of 
the  flat  plate.  Therefore,  great  care  must  be  taken  in  using 
this  type  of  analysis. 

In  order  to  see  what  effect  this  restriction  might  have  on 
the  results  of  the  preceding  section,  the  boundary-layer  theory 
must  be  incorporated  with'  the  heat -conduction  analysis. 

Therefore,  let  us  consider  equations  (9)  and  (12)  and  apply  the 
results  obtained  by  Van  Driest  in  equation  (53).  Assuming  a 
Mach  number  of  unity  equation  (53)  has  been  evaluated  and  is 
presented  in  figure  3  with  Reynolds  number  as  the  parameter. 

Then  by  assuming  a  Reynolds  number  and  free-stream  temperature, 
the  film  coefficient  has  been  evaluated  (refs.  (12)  and  (13)) 
and  is  presented  in  figure  4  for  wall -temperature  ratios  less 
than  unity.  This  is  the  case  of  the  "highly  cooled"  boundary 
layer  and  gives  rise  to  the  most  severe-heating  rates.  The 
film  coefficient  based  on  the  initial  wall-temperature  ratio 
may  then  be  obtained  directly  from  figure  4  and  the  temperature 
distribution  through  the  flat  plate  obtained  by  evaluating  equa¬ 
tions  (9)  and  (12). 

Assuming  the  plate  to  be  initially  at  room  temperature,  the 
temperature  distribution  has  been  calculated  for  various  values 
of  the  Fourier  modulus  and  is  presented  in  figure  5  as  the  family 
of  curves  represented  by  broken  lines.  It  can  be  seen  that 
significant  changes  in  the  wall  temperature  and  therefore  film 
coefficient  occur  even  for  relatively  small  values  of  the 
Fourier  modulus.  In  an  effort  to  see  what  effect  these  changes 
might  have,  the  boundary  layer  was  assumed  to  react  instanta¬ 
neously  to  these  changes  in  wall  temperature.  It  has  been 
reported  (ref.  (14))  that  this  assumption  will  give  satisfactory 
results.  Equations  (9)  and  (12)  have  been  re-evaluated  by 
allowing  the  film  coefficient  to  change  as  the  surface 
temperature  of  the  plate  increases.  This  scheme  produces  the 
temperature  distribution  shown  in  figure  5  as  the  family  of 
curves  represented  by  solid  lines.  Repeating  this  procedure 
does  little  to  increase  the  accuracy  of  the  calculation  as  can 
be  seen  in  figure  6. 


DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

In  the  preceding  section  we  have  discussed  the  effect  of  a 
variable  surface  temperature  on  the  temperature  distribution  in 
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the  one -dimensional  case.  When  considering  the  composite  flat 
plate  it  is  possible  not  only  for  the  surface  temperatures  of 
the  two  materials  to  change  at  different  rates  due  to  the 
difference  in  thermal  properties,  but  the  character  of  the 
boundary  layer  itself  may  change.  Regardless  of  the  cause  of 
the  change  in  film  coefficient  of  heat  transfer,  the  effect  in 
the  two-dimensional  case  may  be  found  in  the  same  manner  as  in 
the  one-dimensional  case.  It  can  be  seen  from  equations  (31) 
that  a  change  in  "h"  has  the  same  effect  as  a  change  in  the 
thermal  conductivity  which  alters  the  Biot  modulus  and  thereby 
changes  the  roots  of  the  transcendental  equation.  It  is 
possible  to  apply  the  iterative  scheme  previously  described  to 
the  two-dimensional  solution  to  determine  this  effect;  however, 
a  simplified  one-dimensional  analysis  based  on  the  maximum  Biot 
modulus  is  sufficient  to  determine  the  maximum  error. 

In  conclusion  it  may  be  said  that  for  the  one-dimensional 
case  considered,  a  maximum  error  in  temperature  ratio  of  less 
than  8  percent  may  be  expected  for  Fourier  moduli  less  than  0.1 
and  less  than  6  percent  for  a  modulus  of  0.03,  by  assuming  the 
film  coefficient  of  heat  transfer  constant  and  equal  to  its 
initial  value.  This  error  is  apparent  only  at  the  surface  of 
the  plate  and  the  temperature  distribution  is  essentially 
unchanged.  The  various  parameters  assumed  in  this  report 
correspond  to  typical  flow  conditions  attainable  in  a 
reflected  shocktube  and  therefore  represent  extremely  severe- 
heating  conditions.  Assuming  the  physical  dimensions  and 
properties  of  a  typical  calorimeter  gage,  the  Fourier  modulus 
of  0.03  would  correspond  to  a  time  of  approximately  300 
microseconds  which  is  well  within  the  range  of  current 
measuring  instruments.  In  wind-tunnel  tests  much  lower  free- 
stream  temperatures  are  present  and,  therefore,  it  is  conceivable 
that  the  constant  film-coefficient  condition  may  be  used  for 
even  larger  values  of  the  Fourier  modulus.  In  reality,  the 
calculated  temperatures  will  always  be  slightly  greater  than 
the  actual  temperatures  and  will  therefore  lead  to  a  conservative 
design  analysis. 

The  equations  derived  in  this  report  for  the  temperature 
distribution  in  a  flat  plate  in  both  one  and  two  dimensions 
have  employed  a  Fourier  analysis  to  determine  the  coefficients 
of  an  infinite  series.  For  small  values  of  the  Fourier 
modulus,  these  series  tend  to  converge  quite  slowly,  making 
hand  calculations  very  impractical.  However,  the  solutions 
presented  here  could  be  incorporated  into  a  machine  program  for 
the  reduction  of  experimental  data  or  prediction  of  heat- 
transfer  parameters. 
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APPENDIX  A 


AN  ORTHOGONALITY  RELATION 

The  decisive  property  of  a  set  of  functions  which  makes 
it  possible  to  determine  one  by  one  the  coefficients  in  the 
assumed  expansion  is  that  the  integral  of  the  product  of  any 
two  distinct  members  of  the  set  taken  over  the  appropriate 
interval  is  zero.  The  set  of  functions  is  then  said  to  be 
orthogonal  (ref.  (3)).  It  may  be  seen  that  the  functions, 

described  by  equations  (35)  ,  are  not  themselves  ortho¬ 
gonal  . 


9.  Jv)  =  cos  -5  costfm  j- 

i 

'  Bp**  COS  C.OS 


(A-la) 

(A-lb) 


Therefore,  in  order  to  perform  the  Indicated  Fourier  analysis 
we  must  obtain  the  orthogonal  functions  associated 

with  the  composite  flat-plate  system.  Let  us  assume  that  the 
set  of  functions  (^W-Cx,y)  are  orthogonal  with  respect  to  some 
"weighting"  function  and  let  us  assume  further  that  function 
to  be  a  constant: 


(A-2a) 

(&■.«“  Ct 

(A -2b ) 

To  prove  the  orthogonality  of  the  <pmat  's  : 

i  a  f  b 

if  ft  ft 

dS  s  +  ))  4.,  J*  Jy  "  0 

(A-3) 

,  5  00  *• 

n  pma  jz  pq  ,  , 

Since  the  satisfy  equations  (26),  so  must  the  <PmtJ  s 

and  ,  therefore: 
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h'' 


(A-4b) 


Now  let  us  multiply  equation  (A-4a)  by  <f>ptf  and  equation  (A-4b) 
by  6  and  subtract  to  obtain: 

( xL O  =*  (A-5) 

Integrating  equation  (A-5)  we  have,  after  expanding  the  right 
hand  side: 

(C-Ofto,,  is  +  ^  - 

5  00  ^ 

j'jc/xo/y  +  JJxt( 

00  ^ 


K 


,  V^^i)  /  /  .((..( 1  /  $  ^ xb±ii 


n  l  ^  x '  ^  y 1 

/VX 

*(  >x 


~  ,  Vh~)) 
V7J 


A  4 


>yV 

(A-6) 


Considering  Just  the  first  Integration  on  the  right  hand  side 
of  equation  (A-6)  we  substitute  from  equations  (A-2)  to  obtain: 

j|  A 


partial  Integration  yields: 


(A-7) 


X  r  CL 


X-  •  o 

Y;i  a 


-4-,  [*..W 

Recalling  our  boundary  conditions: 
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«/ 


<*3., 


‘y  V*. 


=  0 


Y'o 


irL = *.  «-w) 
^riv  =-K(o 


_  K 


(A-10) 

(A-lla) 

(A-llb) 


Therefore,  the  only  contribution  to  the  integral  is  at  x  =  a: 

4 

77^1  .(A  (A-12a) 

o 


Likewise,  from  the  second  integral  on  the  right  hand  side  of 
(A-6)  we  may  obtain: 

*> "ttL-  -^f L i ^  <A-12b) 


At  the  interface: 


3.-,  =  **-&*»,- 

,  and  (A-13) 

"frL=  1'rL 

(A-14) 

Also 

and  (A-13a) 

K'’v  In,,.- ^1... 

(A-I4a) 

Now  substituting  equations  (A-13)  into  (A-12a)  we  obtain: 

k»  jf^i  (A.15) 

0  ••  •  -  —  *■  J«  Ui 


Now  adding  (A-15)  to  (A-12b)  we  obtain: 

*  -If., 


-rr’L;  W*  ^rlj^-x  c,f)jr  <*-*> 
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Since  we  would  like  this  integral  to  be  identically  zero 
everywhere,  let  us  assume: 


(A-17) 


(A-l8a) 

(A-l8b) 


Which  are  the  same  functions  obtained  by  Mayer  (ref.  (4))  for 
the  one-dimensional  composite  slab  problem.  Substituting  the 
appropriate  values  for  Cj  and  C2  into  equation  (A-l6)  we  see 
that  this  integral  is  everywhere  zero.  Therefore,  from 
equation  (A-6)  we  have: 

\\  JS  =  0  (*-i9) 

t 

But  ^  mi  V-  k  mi*  ^ 


Therefore  equation  (A-3)  has  been  satisfied  and  the 
orthogonality  of  the  functions  described  by  equations  (A-2)  has 
been  established.  In  order  to  determine  the  coefficients  of 
the  expansion  lpmm,  we  should  multiply  equation  (A-l)  by  C 
and  integrate  over  the  region  of  interest  to  obtain: 


|cU,Js 

s  s 


( 0 
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'•'»»  *  PI 
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FIGURE  I.  BOUNDARY  LAYER  OVER  A  COMPOSITE  FLAT  PLATE. 
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FIGURE  2.  DETERMINATION  OF  EIGENVALUES 
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FIGURE  3.  LOCAL  TURBULENT  FRICTION  COEFFICIENT  AS 
A  FUNCTION  OF  WALL  TEMPERATURE  AT  MACH  I 
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FIGURE  4  LOCAL  TURBULENT  FILM 

COEFFICIENT  AS  A  FUNCTION 
OF  WALL  TEMPERATURE  AT 
MACH  1 
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FIGURE  5.  TEMPERATURE  DISTRIBUTION  THROUGH  A  FLAT  PLATE 

(M=l,  T  oo=I0.700°R,  R  co=l6.4  X  10® ) 
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FIGURE  6  TRANSIENT  TURBULENT  FILM  COEFFICIENT 
(M  =  I ,  Too  =  l0,700oR,  Roo  =16.4  X  10s) 
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P.  0.  Box  748 
Port  Worth  1,  Texas 

Attn:  Technical  Library 

Gruman  Aircraft  Engr.  Corp. 

Bethpage,  Long  Island,  N.  Y. 

Attn:  Library 

Republic  Aviation  Corp. 

Formingdale,  New  York 
Attn:  Library 

Mart  in -Marietta  Co. 

Baltimore  3>  Maryland 
Attn:  Library 

North  American  Aviation,  Inc. 

Space  &  Information  Systems  Div. 

12214  Lakewood  Blvd. 

Downey,  California 
Attn:  Library 

McDonnell  Aircraft  Corp. 
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